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INTRODUCTION 


This  trajectory  program  is  intended  for  use  when  the  vehicle  can 
be  considered  as  a  point  mass  under  the  influence  of  gravity,  atmospheric 
drag  force,  Coriolis  force  and  centrifugal  force.  The  input  parameters 
have  been  chosen  to  be  as  few  in  number  as  possible  and,  at  the  same  time, 
the  ones  most  often  used.  For  instance,  the  initial  position  of  the  vehicle 
is  determined  by  its  longitude  as  measured  from  Greenwich,  its  latitude 
above  or  below  the  equator  and  its  height  above  the  earth's  surface.  The 
position  of  the  vehicle  in  a  radar  system  is  calculated  from  a  knowledge  of  the 
longitude  and  latitude  of  the  radar  site. 

The  equations  of  motion,  coordinate  transformations  and  auxiliary 
computations  are  contained  in  the  main  body  of  the  report.  Appendix  A 
contains  the  numerical  integration  procedure  used  to  solve  the  equations 
of  motion.  Appendix  B  contains  a  list  of  the  Fortran  symbols  and  corres¬ 
ponding  mathematical  symbols,  a  complete  listing  of  the  Fortran  state¬ 
ments  and  a  flow  diagram  of  the  program.  The  input-output  quantities 
are  also  defined  in  Appendix  B. 

The  author  would  like  to  thank  Mr.  Thomas  J.  Kroupa  III  for  his 
assistance  in  programming  the  equations  for  an  IBM  1410  computer. 
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EQUATIONS  OK  MOTION 


The  trajectory  of  the  vehicle  is  referenced  to  a  r  1  ght  -  handed 
rectangular  co-ordinate  system,  x,  y,  z,  rigidly  connected  to  the 
rotating  earth  and  with  the  origin  at  the  earth's  center.  The  z-axis 
is  along  the  earth's  polar  axis  and  the  xy  plane  is  in  the  plane  of  the 
equator  with  the  x-axis  located  at  the  meridian  of  Greenwich.  (See 
Figure  l ) . 

The  forces  acting  on  the  vehicle  in  this  system  are  gravitational, 
air  resistance,  Coriolis,  and  centrifugal.  With  the  assumption  that  the 
force  due  to  air  resistance  varies  as  -k  ,  the  equations  of  motion 
along  each  ot  the  co-ordinate  axes  are:  (Reference  1) 


- Zl — - -  _  K  k(x?  f  +  v.1 )  -  l  £  y^>  t  u/  \  (I) 

(xz  ♦  y"  ♦  /M  ' 


-K  v 

i 

-  k  y  ( xz  *  yz  *  i  : 

y  '  (x4  «  V-  ♦ 

...  -  “P"  * 

1 

-  k  z(xz  ♦  yz  *  >*)' 

(xi  +  yz  f 

■/}  w 

.  x 

where  x,  y,  z  -  — y- 

dr 

nr 

respectively,  and 
dT 

d  x 

x-  y’  7  -  hT  ' 

±L 

dt 

respectively 

dt 

Polar  Axis 


Figure  1 

Relationship  Between  the  Earth-Fixed  x,  y  ,  z  System 
and  the  Radar  Xj  ,  yj  ,  Zj  System 


5 


^  =  1*0  Rc' 

g^  =  acceleration  of  gravity  at  sea-level 
Re  =  mean  radius  of  the  earth 

k  =  1  K0  p/P 

p  =  atmospheric  density 

p  is  computed  as  a  function  of  altitude  by  a  subroutine 
based  on  the  ARDC  Model  Atmosphere,  1959  (Ref.  4) 

ft  =  W/CqA  ,  the  ballistic  coefficient 

W  =  the  weight  of  the  vehicle 

Cp  =  drag  coefficient 

A  =  reference  area 

u>  =  earth's  rotation  rate 

Equations  (1),  (2),  and  (3)  were  numerically  integrated  by  the 
fourth  order  method  of  Runge-Kutta  as  outlined  in  Appendix  A.  To 
start  the  integration  procedure,  a  point  in  the  7 -dimens ional  configuration 
space  t  x  y  z  k  y  z  must  be  known.  This  point  is  determined  from  the 
usual  earth  referenced  trajectory  parameters  (speed,  altitude,  latitude, 
longitude,  bearing  angle  and  re-entry  angle)  by  the  following  transforma¬ 
tions.  (See  Figure  2) 

x  -  [  Re(4>)  +  H]  cos  4>cos  0 

y  =  I  Re W>)  +  H]  cos  ^sin  0  (4) 

z.  =  l  Re (4>)  +  H]  sin  0 


4 


V 


(5) 


* 

X 

y 


-ain  0 
cos  0 
0 


-sin  <}>  cos  0 
-sin  4>  sin  0 
cos  4> 


sin  4> 


0 

Vx 

xm 

e 

V 

Vm 

Vz 

zm 

where  Vv 

xm 


0 

4> 

6 

4* 


V  cos  6  sin  ^ 

V  cos  6  cos  ^ 

V  sin  ^ 

longitude  of  the  vehicle 
latitude  of  the  vehicle 

re-entry  angle  (positive  upward  -  see  Figure  2) 

bearing  angle  (positive  clockwise  from  north  -  see 
Figure  2) 


VX  •  Vy  - 

xm  y  m 

Vz^  =  the  components  of  V  in  the  vehicle's  local  reference 
system  -  see  Figure  2 

H  =  altitude  above  the  earth's  surface 

Reto>)  =  radius  of  the  earth  at  latitude  4> 

i 

Re(4>)  =  20855967(1  -  .00672267  cos2  <$>)*  2 


Equations  (4)  and  (5)  are  integral  parts  of  the  computer  program 
The  vehicle's  position  at  any  time  during  the  integration  is 
given  by: 

i 

H  =  (x*  +  yl  +  z*)2  -  Re(4>) 


(6) 


6 


z 


(?) 


4>  =  tan  “  - 

(x*  + 

where  if  z  is  positive, 
negative,  -90°  4>  < 


0  <  4>  <  90°  (north  latitude),  and  if  z  is 

0  (south  latitude)  . 


0 


0  =  tan'1  1  (8) 

X 

To  remove  the  ambiguity  from  G,  east  longitudes  were  chosen  to  be 
positive,  and  west  longitudes  negative.  The  value  of  9  is  determined 
from  Subroutine  QUAD  by  the  following  scheme: 
y  x  0 

+  +  O*<0<9O° 

+  -  90°  <  0  <  180° 

+  -90*  <  0  <  0° 

-180°  <  0  <  -90° 


Ground  range  (S)  is  defined  to  be  the  distance  traveled  from 
the  initial  point  along  the  earth’s  surface  and  is  computed  in  increments 
as  follows: 

Let  R  =  i  x  +  j  y  +  k  z  at  time  t 

and  Rq  =  i  Xq  +  j  yQ  +  k  zQ  at  time  t  -  At 

Then: 

|RqxR|  =  |R0|  |  R  |  sine 

where  4  1®  the  angle  between  the  two  vectors.  Since  the  computation 
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interval  At  is  very  small,  R  will  differ  very  little  from  RQt  and 
and  sin  §  = 

Thus , 

|  =  |Rq  x  R|/|Rj  |R|  .  and 


AS  =  Re  4 

AS  is  summed  at  the  end  of  each  computation  interval  to  give  S. 
T*  e  re-entry  angle,  6  ,  is  given  at  any  time  by: 


(9) 


R  •  V  =  |  R|  |  V|  cos  (tt/2  -  6) 

=  1  R  I  I  Vi  sin  6 


6  =  sin'1  R  •  V  /  I  R|  I V 


(10) 


where: 


and, 


R  =  ix+jy+kz 


V  =  ix+jy+kz 


6  is  defined  to  be  positive  when  above  the  local  horizontal  and  negative 
when  below.  (See  Figure  2). 
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TRAJECTORY  PARAMETERS  IN  A  RADAR  REFERENCE  SYSTEM 


The  co-ordinate  system  xj  ,  y{  ,  with  origin  O,  at  the  radar  is 
defined  as  follows:  (See  Figure  1) 

xi  »  Yi  »  zi  are  co-ordinate  axes  with  origin  O,  at  the  surface  of  the 
earth,  the  xj  ,  yi  plane  is  perpendicular  to  a  radius  vector  drawn  from  the 
center  of  the  earth  and  the  Zi  axis  is  along  the  radius  vector,  the  positive 
direction  or  xj  and  yj  are  taken  to  be  due  east  and  due  north  respectively. 

For  a  station  at  latitude  4>r  an(*  longitude  0r,  the  co-ordinates  of 
Ol  in  the  earth-fixed  x,  y,  z  system  are: 

xr  =  Rer  cos  c 08  ®r 

yr  =  Rer  cos  «|)r  sin  0r  (11) 

zr  =  Rer  sin  4>r 

where  Rer  is  the  value  of  Re(<t>)  at 

Using  the  standard  equation  for  translation  and  rotation  of 
co-ordinate  axes,  the  following  relationship  between  the  two  systems 
is  obtained.  From  x,  y,  z  to  xj  ,yj  ,  Zj  : 

t 

t 

X]  -sin  0r  cos  Br  0  I  x  -  xr 

yj  =  -sin  cos  0r  -sin  $r  sin  cos  ♦r  Y  *  Yr  (12) 

|  Z|  cos  4>r  c o 8  0f  cos  4>r  sin  0r  sin  4>r  z  -  zr 

i 

i 

i 


i 
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i' rum  x  i  ,  y  [  ,  z  j  to  x,  y  ,  z: 


X 

-sin  0 

-sin  4>r 

cos  0r 

y 

- 

cos  0|. 

-sin  4>r 

sin  0r 

7. 

1 

! 

0 

COS  < 

cos  4>  cos  0r 

i  X! 

xr 

cos  bin  0r 

t 

yi 

4- 

(13) 

sin  <\)  . 

1 

! 

*r 

hi  tlu*  computer  program*  these  tr anstor MiaUons  are  executed 
I) y  the  subroutines  GOOD  and  GOOD!.  Thus,  if  one  wishes  to  define  the 
radar  system  in  some  other  manner,  only  the  subroutines  will  have  to 
be  ch.tngi  G 

After  the  vehicle’s  position  has  been  transformed  from  the 
x,  y,  z  system  to  the  x  lt  yj  ,  Z\  system,  the  slant  range,  azimuth 
angle,  and  elevation  angle  are  computed  as  follows: 


_i 


Hi 

=  (x,* 

+  y/  + 

(14) 

—  tan 

[  Zi  /(X I1  +  y/)*  ! 

(15) 

A* 

=  tan'1 

(- ) 

;i6? 

'  yi  / 

E f  ranges  from  0#  to  90°  and  is  positive  if  the  vehicle  is 
above  the  horizon. 

Az  ranges  from  0*  to  360°  and  is  measured  positive  clockwise 
t  r i >rn  north . 

R .  ,  E|,  and  Az  are  computed  in  subroutine  RAE.  The  comments 
made  about  GOOD  and  COODI  also  apply  to  RAE. 
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Hi  can  be  expressed  in  the  x,  y,  z  system  as: 


Hi  =  i  (x  -  xr)  f  j  (y  -  yr)  f  k  (/.  -  y.r) 


and  the  velocity  vector  in  the  name  system  is: 


V  =  i  x  +  j  y  +  k  7, 


Using  these  two  equations,  the  velocity  aspect  angle  and  the  range  rate 
can  be  computed. 

The  velocity  aspect  angle  (angle  between  the*  radar  line  of 
sight  and  the  velocity  vector;  it  given  by: 


V 


cos 


R 1  .  V 

Ri  l  | V  ' 


17) 


The  component  of  |V|  along  R1  is  the  range-rate  (Rj).  Thus, 

R|  ;  |  V  |  cos  v  (18) 

Since  R|  is  negative  when  the  vehicle  is  approaching  the  radar, 

Y  is  chosen  to  range  from  0°  to  180°.  y  =  ®  °  when  the  vehicle  is 
going  directly  away  from  the  radar,  and  y  =  180°  when  *he  vehicle 


1  i 


is  headed  straight  in. 


CONCLUSIONS 


The  output  of  the  computer  program  has  been  compared  to  actual 
radar  data  and  found  to  be  in  good  agreement.  It  is  felt  that  the  program 
will  be  useful  for  generating  theoretical  slowdown  curves  and  for 
determining  range,  range  rates,  and  look-angles  from  arbitrary  radar 
locations . 
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APPENDIX  A 


NUMERICAL  INTEGRATION  OF  THE  EQUATIONS  OF  MOTION 

The  method  described  below  in  the  classic  fourth  order  procedure 
of  Runge-Kutta  (Reference  2).  Only  one  point  on  the  integral  curves  is 
needed  to  start  the  integration,  and  with  the  aid  of  high  speed  computing 
machines,  any  degree  of  accuracy  can  be  achieved  by  choosing  a 
sufficiently  small  increment  of  the  independent  variable. 

Writing  equations  (1),  (2),  and  (3)  as: 


x  =  fj  (t,  x,  y,  z,  x.  y,  z) 

(1A) 

y  =  f2  (t,  x.  y,  z,  x.  y,  z) 

(2A) 

z  =  f3  (t,  x,  y,  z,  x,  y,  z) 

(3  A) 

The  integration  proceeds  as  follows. 

Let  t  take  on  an  increment  At;  then  x,  y,  z,  x,  y,  and  z  receive 

increments  Kj  ,  Kt ,  K3  ,  K4  .  Ks  ,  and  K6  respectively. 

Ki  =  1/6  (k„  +  2k2,  + 

2k3i  +  k41) 

(4A) 

K2  =  1/6  (ku  +  2kw  + 

2k3Z  +  k42) 

(5A) 

K,  =  1/6  (k,,  +  tk2s  + 

2  k33  +  k43 ) 

(6A) 

I  3 


K*  =  1/6  (k,4  +  2k24  +  2  k34  +  k44)  (7  A) 

K5  =1/6  (k i5  +  2  k25  +  2  k35  +  k45 )  (8A) 

K6  =  1/6  (ki6  +  2k26  +  2  k36  +  k46)  (9A) 

k  |  j  =  xAt  (16A) 

k j2  =  y  At  (1  1  A) 

ku  =  zAt  (12A) 

kj4  =  fj  (t,  x,  y,  z,  x,  y,  i)  At  (13A) 

k»  =  f  2(t,  x,  y ,  z,  x,  y,  z)  At  (14A) 

k36  =  fj  (t,  x,  y,  z,  x  .  y,  z)  At  (1  5A) 

k23  =  (x  +  |k|4  )  At  (16A) 

k22  =  (y  +  2k,5  )At  (17A) 

k21  =  (z  +  2Lku)At  (18A) 

k^  =  f|  (t  +  i  At,  x  +  |kn  ,  y  +  ik22  ,  z  +  2kjj  ,  x  +  • 

y  +  |ki4  ,  i  +  ik16)  At  (1  9 A) 
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(20  A) 


I 

I 

I 

I 

I 

l 

I 

I 

I 

I 

\ 

I 

I 

l 


k25  ^2  (t  ^  2  At,  X  +  2k  jj  ,  y  +  2ki2  >  z  2k?.3  ,  X  +  2k14  • 

y  +  |k15  ,  z  +  jkl6  , )  At 


k2fe  ~  ^3  ^  2  At,  x  2k  jj  ,  y  +  ^  kj2  >  z  2kj3  ,  X  +  2kj4  , 

y  +  xki5  ,  z  +  |k16,)  At  (2  i  A) 

k3i  =  (k  +  |k24)  At  (22A) 

k32  =  (y  +  ik25)  At  (23A) 

k3j  =  (z  +  jk26  )  At  (24A) 

kn  =  ft  (t  +  j  At,  x  +  £k2,  ,  y  +  2  k22  ,  z  +  Jk23  .  x  +  |k24  . 

V  +  £k25  ,  z  +  £k26)  At  (2 5 A) 

kj5  =  h  (t  +  2  At.  X  +  2  k21  .  y  +  2k22  .  z  +  2  k23  .  x  +  Ik24  , 

y  +  2k25  ,  z  +  ]k26)  At  (2 6 A' 

k15  ~  ^3  M  i  At,  x  ^  2 k  21  I  y  +  2 k22  »  z  ^  2k23  •  x  ^  lk24  > 

y  +  ik25  ,  z  +  2k26)  At  (27 A) 

k4i  =  (x  +  kJ4)  At  (28A) 

k42  =  (V  +  k35  )  At  (29A) 

k4,  =  (z  +  k,t)  At  (30A) 

k44  -  m  (t  +At,  x  +  k}(  ,  y  +  k3  2,  z  +  k3i  ,  x  +  kJ4  ,  y  +  k35  , 

z  ♦  k36  )  At  (31  A) 


IS 


y  +  k,5  , 


k4S  =  I',,  (t  t  At,  x  +  ku  ,  y  +  ku  ,  z  +  k ,  x  +  k34  , 

z  +  kj6)  At  (32 A) 

k4b  -  t,  (t  »  At.  x  +  k,,  ,  y  +  kj2  ,  z  +  k)3  ,  k  +  kJ4  ,  y  +  k>s 

z  3  kJ6)  At  (33 A) 

On  the  first  pass  through  equations  10A  to  3  3A,  the 
variables  t,  x,  y,  z,  x,  y,  z  will  have  their  initial  values.  After 
equation  3  3 A  has  been  executed,  t  is  incremented  by  At,  x  by  Kj  , 
y  by  Kz,  z  by  K3 ,  x  by  Kt,  y  by  K5 ,  and  z  by  Kb  and  the  procedure 
beginning  at  10A  is  repeated. 

During  portions  of  the  trajectory  where  the  acceleration 
is  small  At  may  be  chosen  fairly  large  (around  \  sec),  but  when 
the  acceleration  is  large,  At  must  be  small  (around  1/100  sec). 

Equations  fj  ,  il  ,  and  f*  are  evaluated  in  the  program  by 
subroutine  FU12  3. 
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List  of  FORTRAN  Symbols 
Inputs  to  the  Program 
Output  of  the  Program 
Sample  Print-Out 

List  of  FORTRAN  Statements 

Flow  Diagram 
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LIST  OF  FORTRAN  SYMBOLS 


FORTRAN  Symbols 

Mathematical  Symbols 

AKl ,  AK2,  AK3 

K,  ,  K2l 

AK4,  AK5,  AK6 

k4.  k5> 

k6 

AKl  1 ,  AKl  1,  AK  13 

kn  ,  kJ2 , 

-  ki3 

AKl 4,  AKl  5,  AKl 6 

ku  .  kls 

-  kn 

AK21,  AK22,  AK23 

k21  ,  k22 

►  k23 

AK24,  AK25,  AK26 

kz4  .  k2s 

■  k2(> 

AK31 ,  AK32,  AK33 

kj>  .  k}2 

,  k33 

AK34,  AK35,  AK36 

k}4  ,  k„ 

■  kl6 

AK4 1 ,  AK42,  AK43 

(M 

.  k43 

AK44 ,  AK45,  AK46 

k*4  .  k45 

>  k4t 

ALA 

*r 

ALO 

0 

r 

AT 

X  2  +  y 

2  +  z  2 

AZ 

A 

z 

BC 

p 

BETA 

4- 

DEL 

* 

ip  t 

r#  Ij 

I)T 

At 

Gamma 

V 
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FORTRAN  Symbols 


Mathematical  Symbols 


H 

PHI 

RAN 

RER 

RO 

RV 

R 1 

RR1 

T 

THETA 

V 

VT 

VX,  VY,  VZ 

X.  Y,  Z 

XD,  YD.  ZD 

XDI),  YDD,  7.DD 
XR,  YR.  ZR 

XI,  Yl,  7.1 


H 

4> 

S 

R  (4>)  at  4>  =  <t> 
e  r 

P 

H  +  R  (<t>) 
e 

Ri 

Ri 

t 

o 

v 

(x2  +  y 2  +  z2)2 


X.  Y.  Z 

•  9  I 

X.  Y,  Z 


X,  Y.  Z 

X  ,  Y  .  Z 
r  r  r 

X,.  Y,,  Z: 
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INPUTS  TO  THE  PROGRAM 


The  initial  conditions  are  read  into  the  computer  from  four 
input  cards  containing  the  following  information: 

Card  #1  BC,  DEL,  BETA,  V 

BC  =  the  ballistic  coefficient  (W/CqA)  in  Ibs/ft^ 

DEL  =  the  re-entry  angle  (negative  when  re-entering)  in  degrees 

BETA  =  the  velocity  bearing  angle  (positive  clockwise  from 
north)  in  degrees 

V  =  magnitude  of  the  velocity  in  ft/sec 


Card  #2  THETA,  PHI,  H 

THETA  =  longitude  of  the  vehicle  in  degrees.  If  the 
longitude  is  given  as  0  degrees  west  of 
Greenwich,  change  to  360°  -  0. 

PHI  =  latitude  of  vehicle  in  degrees  -  input  as  positive 
when  above  the  Equator  and  negative  when  below 

H  =  altitude  above  the  earth's  surface  in  feet 


Card  H  ALO,  ALA 

ALO  =  longitude  of  radar  site  (input  in  the  same  manner 
as  THETA) 

ALA  =  latitude  of  the  radar  site  (input  in  the  same  manner 
as  PHI) 
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Card  #4 


DT,  HEND,  N 


DT  =  increment  of  the  independent  variable  time 

HEND  =  altitude  at  which  it  is  desired  that  the  compilation 
be  halted 

N  =  print  rate  -  controls  the  number  of  times  through  the 
integration  loop  before  printing.  For  instance,  if  DT 
is  1/10  sec  and  trajectory  data  is  desired  at  1  sec 
intervals,  then  N  should  be  read  in  as  10. 
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OUTPUT  OF  THE  PROGRAM 


The  quantities  shown  on  the  sample  print-out  are  defined  as 


follows : 


TIME  =  Elapsed  time  in  seconds  from  initial  point 

LATITUDE  =  Latitude  of  the  vehicle  in  degrees  (positive  when  north 
of  the  equator) 

LONGITUDE  =  Longitude  of  the  vehicle  in  degrees  (positive  when 
east  of  Greenwich) 

ALTITUDE  =  Height  above  the  earth's  surface  in  feet 

TOTAL  ACCELERATION  =  Absolute  value  of  the  acceleration  (ft/sec2) 

in  the  earth-fixed  reference  system 

TOTAL  VELOCITY  =  Absolute  value  of  the  velocity  (ft/sec)  in  the 

earth-fixed  reference  system 

GROUND  RANGE  =  Distance  traveled  over  the  earth's  surface  from 

the  initial  point  (nautical  miles) 

RE-ENTRY  ANGLE  =  Angle  between  the  velocity  vector  and  the  local 

horizontal  (degrees) 

X  RADAR 

Y  RADAR  =  Co-ordinates  of  the  vehicle  in  the  radar  reference 
Z  RADAR  system  (feet) 

AZIMUTH  ANGLE 

ELEVATION  ANGLE  =  Radar  look-angles  in  degrees  (azimuth  is 

no.  of  degrees  clockwise  from  north,  and 
elevation  is  no.  of  degrees  above  the  horizon). 

GAMMA  =  Angle  between  the  radar  line-of-sight  and  the  velocity 
vector  (degrees) 

SLANT  RANGE  =  Distai.*  c  from  the  radar  to  the  vehicle  (feet) 

RANGE  RATE  =  Rat r  af  change  of  the  slant  range  (ft/sec) 
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LIST  OF  FORTRAN  PROGRAM 


PAIN  PROGRAM 


FORTRAN  RON 
BCP 

LINE  *  C 
IPAGE  *  1 

RE  AC  100, BC,OEL, BETA,  V 
READ  100, THETA, PHI ,H 
RE AC  10C,AL0,ALA 
READ  101 ,DT ,HENO,N 
PRINT  1C7 
PRINT  IC6 

PRINT  1 C2 ,H , V  , DEL. BC, BETA, PHI , THETA 
RAN«0. 

I  (HO. 

MOO 

DEL  #DEL  •.01  745 
BE TA#BETA*. 01745 
THETA#THETA*. 01745 
PH  I #PH I  * . 0 1 745 
AL0l»AL0«. 01745 
ALA#ALA«. 01745 
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LIST  OF  FORTRAN  PROGRAM 


MAIN  PROGRAM 


RV#h&208‘i‘S9ft7./SQRTFSl.-.006722fi7*C0SFl<PHIn**2n 

KERrf208,;‘j<>fc7./SQRTFSl.-.006  7226  7*CGSF*ALAa**2n 

XR#RER*COSF¥ALAa*COSF*ALOn 

YR(KRER*COSFHALAa*SINFXALOa 

ZR#RER*SINF*ALAa 

VX#V*COSF«CFLD*SINr*BETAa 

VY#V«COSF*CELQ*COSF*BET Ad 

VZ#V»S  INFtCELn 

CALL  r.ncn*xn  ,YO  ,ZO  ,VX,VY,VZtTHE  TA,PHIa 
XKRV*CO«;F%PH|o»COSF«THE  tad 
Y«RV*CCSFXPHIa*SINFlTHE  TAD 
Z  HR V*S I NFtPH I o 
xo*x 

YO»Y 

10*1 

RVOiRV 

CALL  ALT*H,RO,A,BtC.D,E,F,Gc 
RC0RGM2.17A 
AK 1 1 #XO»DT 
AK  1  2(KYD*DT 
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LIST  OF  FORTRAN  PROGRAM 


MAIN  PROGRAM 


AK  1 

3  A  Z  0  •  D  T 

CAL 

L  F  U I 23% AK  1  4 

9  AK 1 5 i 

,AKl<S, 

,X,Y,Z  tXD, YD»  ZD,RO,BC ,DTa 

AK  2 

1 *%XD£. 

5  •  AK  I 

An«D  T 

AK  2 

2*SYC£. 

5  •  AK  1 

5n#D  f 

AK  2 

3#%ZD£. 

5  •  AK  l 

6a#0T 

OCAL 

L  FU123IAK24 

•  AK25, 

i  AK26 i 

,X£.S«AKI  l ,  Y£.r»*AKl2,  Z£.5*AK13,  XO£.r»»AK  14 

I  ,Y0£.*i*AKI5,Z06.5*AKI6fR0,8C,DTo 
AK31*%XC£.‘j*AK24d«DT 
AK3  2H»1:YD£.,5»AK2‘5d*DT 
AK33mZD£.«S*AK2ftn*DT 

OCALL  FU123*AK34,AK35,AK  36 ,  X£ . •»•  AK2 1 ,  Y£ . S *AK22,  Z£ .5* AK2 3,  XD£.*»*AK24 
1 , YD£.5»AK25, ZD£.5*AK2A*Rn,0C,0To 
AK4lf»%XD£AK34n*DT 
AK42(TtYC£AK3So*DT 
AK43A*ZC£AK36n*OT 

OCALL  FU123IAK44»AK45*AK46,X£AK31tY£AK32»Z£AK33#X0CAK34»YD£AK3^»Z0£ 

1AK36,RO,BC,DTO 

AKlftSAKl  !£2.*AK21£2.*AK31£AK4lo/f . 

AK2(**AK12£2.«AK22£2.  *AK  32CAK42C 
AK  3 0 YAK  l  3£2.»AK2  3£2.*AK33£AK4  3ti/4. 
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LIST  or  FORTRAN  PROGRAM 


MAIN  PROGRAM 


AK  4  Af  *  AK  \UL2  .  •  AK  24  £?  .  *  AK  34CAK44Q/6. 

AKSurtAK  l^C?.  »AK?*iG2  .*  AK  3^CAKAr)Q/(S. 

AK6#*AK  INC 2.  »AK?Gt2  .  *AK  3NCAK4M3/N. 

XAXCAK1 

YMYCAK? 

ZHZGAK3 

XCMXCCAK4 

YCHYOtAKS 

ZCUZCCAKA 

RV*SCRTF*X*X£Y*YGZ«  Zo 

HI*KV-208'>'>9fr7./SORTF%l.  -  .006?226  7**X«XGY*Yn/RV*«2a 

lllAHZ  •YC-Y*ZOn*«2£%X«ZO-Z*XCu»»2£*Y*XO-X*YOa#«2 

C*SURlF»Blo/*Rv»RVOQ 

RAN#C«2.O70SO*>E£O7/6O7A.1£RAN 

TUTtCT 

M«MC1 

IFtN-MD2,2,  l 

2  CALL  OUAntTHE TA»PHI t XtY  ,Zo 
MHO 

VT(TSCRTF1!XC*XOGYO«YR£ZC*ZOq 
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LIST  or  FORTRAN  PROGRAM 


MAIN  PROGRAM 


f.Al  l  T  U  1  2  3‘JXOD,  YOU,  ZOO.  X,  Y,  Z*XO»YO,ZD*RQ»BC» 

Ai«soRTF‘4;xon«xcncYDO*YDno/DO*/r:rn 

l)#*X»Xi:CY»YnCZ»ZOa/HRV«VTn 

OUTSORT!  *1  . -l)*Qo 

OK  <TAT  ANF*C/01a/.Cl  74*>  33 

OAK  COCOUXl  ,Yl  ,Z  l  ,X-XR,  Y-YR,Z-ZR,ALO,ALAa 
CAll  RAf *R 1 ,AZ ,Ft ,X 1 , Yl  ,Z  lu 
RR  1  )***X-XRD*XOt«Y-YRa#YOt*Z-ZRD*ZOa/Rl 
CGAMA#RR1/VT 

sgamaascrtfi:i.-cgama**?q 
If XRR In  12, 12, 13 

I  3  GAMAAATANFtSUAMA/CGAMAD/.Ol  7<.S 

GO  TO  1  *» 

1  2  GAMA#18C.-ABSFtATANF*SGAMA/CGAMAn/.Ol7<,Sa 
I S  CON  1  I NUF 

PRINT  200, T , PH  I , AZ 
PRINT  201, H, THETA, EL 
RRINT  202, AT, DEL, GAMA 
PRINT  11U.X1.M,/' 


PRINT  203.VT.R1 


LIST  OF  FORTRAN  PROGRAM 


PAIN  PROGRAM 


PRINT  2C4.RAN.RR1 
l  INF  I i  l  INt  fc  I 
IFUINE  -  fcol?34,  U23, 1234 
142J  LINE  X  C 

I  PAGE  K  I  PAGE  t  1 
PRIM  1532.1  PAGE 
1ST?  FCRPATXIH1  »  75X.4HPAGE.  I  50 
12T4  CONTINUE 

IFtH-HCNODA.A,  1 
<.  PRINT  r*CO 
STOP 

100  FCRMATX8F 15.8o 

101  FCRKATX2E15.8. 13 

; 0 1  FCRPATSIHl .  5  fx . 1 8H INITIAL  CCNOIIIONSn 

lU^OFCRPATt IHK, 25X, 3HAL T » 5X » 8HVEL0C I T Y » 3X . 10HRE-ENT  ANG.7X.6HBAL  C0.5X 
I .8HBEAR  ANG.9X.4HLAT0. 9X.4HL0N0O 

102  FORPATXIH  » 17X.lOEI3.5o 

2000F0RMATX /// » 18X  »  AH? I  PE  » 1 8X . £ I2.5«6X tSHLAT I TUOE  » 9X »E 12. 5 » 5X » 1 3HAZ I HU 

1 TH  ANGLE. 5X.E12. 50 

20 1 OFCRMATX I H  . 15X.8HALT1 TUDE. 1 AX »E 12. 5 .5X.9HL0NGI TUDE » 8X» E l 2. 5. 5X # I5H 
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LIST  OF  FORTRAN  PROGRAM 


MAIN  PROGRAM 


1ELCVAT  !CN  ANGL E  ,  3 X ,  E  1  2.  r.u 

2020FCRMATS1H  ,  1SX  ,  1  8HTCTAL  ACC  E  IE  R  A  T  I  ON  ,  A  X  ,  E  1  2 . 5 ,  ft  X  ,  1 AHRE-ENTR  Y  ANGLE 
1  ,  3  X,  E  1  2.  •» ,  •>  X,  r.HG  AMM  A,  13  X,  El  2.5a 

2030FURMAU1H  ,  1 5X ,  1  Ah  TOT  AL  VE LCC l T Y , BX , E 1 2. S , 6X , 1 1 HSL ANT  RANGE, 6X,E12 
I.Scj 

20 A  FORMATS  IH  ,  1 5X , 1 2HGR0UND  « ANGE  ,  1  OX , E  1 2 . ** , 6X  ,  10HR ANGE  RATE,  7X  ,E12.*> 
lo 

till  FERMATtlH  ,15X,RHX  RAOAR  ,  I  AX , E  1 2. 5 , 6X , 8HY  RAOAR  ,  9X , E 1 2 . 5 , 5X , 7HZ 
IRADAR, I  IX, E 12.50 

SOO  Ff.RMATXlHl  , f OX,  10HEN0  OF  JOB  /////////o 
END 
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LIST  OF  FORTRAN  PROGRAM 


SUBROUTINE  COOD 


bop  conn 

SUBROUTINE  COOD* A f B ,C , D , E , F , 0 , Po 

A#-G»SlNF*Gn-E*SINF*Pa*COSF*OntF«COSF*Pn«COSF*Oa 

B#D*COSF*Od-C*SINF*Po»S INF*CotF*COSFXPn»SINF*On 

C#t  *COSI %PoEF  «S I NF  XPa 

RETURN 

ENU 
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LIS!  OF  MM!  KAN  PROGR  A  M 


SUBROL T INC 


bc:p  UUAD 

SUHKUJ  I  I  N't  QUA01IA*B#Xf  Yf/n 
A  #  A ! ANF  YY/Xo/.O!  7AS 
U *Yv] . 2, 2 
'  lF«Xn39A9?0 

3  A  H  -IRC.  C  A 
GU  IU  2  C 
<♦  AM90. 

GC  TO  20 
/  U*XO6fAf?0 

s  a  #  i  e  o  •  c  a 

GG  I C  2C 

2  0  B#AT ANFXZ/SSURTFSX«X£Y*Ynno/.01  7 AS 
KI  1  UR  N 

f  NO 


*  } 


(JUAf) 


LIST  OF  FORTRAN  PROGRAM 


SUBROUTINE  RAF 


ROP  RAE 

SUBROUTINE  RAE*A,8»CtXl ,Y1 , L lu 
A0SCRTFSXI«X1  t  Y  l  •  Y  1  L  ZWln 
G*AIANF*Zl/S(jRTn;Xl*Xl  L  Y  1  •  Y  1  00/ .  0  l  74S 
H#ATANF*X1/Yln/.01 74S 
IF«XlQ|,2*2 

1  IF*Y1d3,4,S 
\  ItllRO.EB 

r.f  il  is 
B  #  9  0  • 

G(.  TO  IS 
r>  1!#37>C.£B 
GO  TO  IS 

2  IF»Yloft,A,  IS 
f>  R4180.ER 

GO  IC  IS 
lr>  RETURN 
END 
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LIST  Of-  TOKfKAN  PMOC-RAM 


SUHRCU  T I Nt  conn  I 


HCP  COODI 

SLBRCUT INt  COOO!*A,B,C.O,EtF,n,Pn 
A#-n*S  IM'tCntE-COSFSOa 

B(»-C»SiNF«PD»COSF*OD-£»SINF*Pn«SINF*Oa{.F  •COSFSPa 
C#C*COSF*Pa«COSF*Coee*COSF*Pn*S  INF*m».F«SINF»Pn 
RETURN 
END 
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LIST  QF  FORTRAN  PROGRAM 


SURROL  T INE  FUI23 


BCR  FIJI  23 

SUBROUTINE  FU 1 2 3%A , B *C , X , Y , l , XD , YD , 20 , RO , BC ,  IH D 

C2#  1  .  3899909  IE G  16 

U H.  72918296E-04 

C  3  rfl 6 . 0  8  7 • RC / BC 

OH*X»XGY*YGZ*Zn**l  .  S 

D2#SCRTF%XD*XD£YO«YO£ZD*/Oci 

Ad(J:-C?«X/0l-C3*X0*C)2C2.»O«YC£U»lJ*Xn»DT 

Bl*t-C2»Y/ni-C3*YO»n2-2.  •U»XCCU«U*Yn»OT 

C»X-C2*Z/ni-C3*ZD»02n*DT 

RETURN 

ENU 
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LIST  OF  FORTRAN  PROGRAM 


SUBROUTINE  ALT 


RCP  ALT 

SUBROUT INE  AL I , RO , PR , CF , F P , T E M , GR A , WM , SOSa 
lF*P-4OCOOC.n30, 30, 31 
31  RC#0. 

RETURN 

30  IF*H-40CO0.Dl,2,3 
1  X#H/  lOOCO. 

P3#  -0.S672184CE-C3 
P2#  -0.9S8C8049E-02 
PI#  -0. 3734  7339EC00 
PO#  0. 191 73S74ECC2 
T  3#  0. 12881606ECC1 

T  2#  -0.E0S34827Ef.Cl 
Tl#  -0. 2881 3737ECC2 
10#  O.S17723BSLCC3 
GO  TL  SC 
?  CONTINUE 

3  IFXF-BCC00.a4tS,<S 

4  CONTINUE 

5  XHSF-40000.0/ 1000C. 
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LIST  OF  fCRIRAN  PROGRAM 


SUBROUTINE 


"  1#  0.  T9  193A40E-0S 
IV#  0.2S19732AE-03 
>M  0-0.4  7  8  H  3  A  2  3  E  6  0  C 
<’<)#  0.  1  7A8707.9EG02 
’  T#  0. 

!2#  C. 

!  1  #  (J  . 

TUI  C. 3B999000EE03 
Gl  It  SC 

f  1  I  fH-l<SC000. 07,8,9 
7  XK^H-ROCOO.o/lOOOC. 
I*  T#  - 0.  I  9A231  33E-03 
l'2#  0.99IS76253E-02 

1*1#  -O.A8O90213EE0O 
‘‘0#  0.  t  SS7S633EE02 

M#  -0.  109A28A8EC00 
12#  0. 1231 1 8BAEE01 

>1#  0.  12AKS9I  3EC02 

10#  0. 3892S80SEC03 

Ml  If)  SC 


AL  I 
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LIST  OF  FORTRAN  PROGRAM 


SUBROli  TINE 


8  CON  I  I NUF 

9  IF?H-1 75C0C.U10,  lit  1? 

10  CONTINUE 

11  X**h-l*OOOC.a/lCOCO. 

P  3#  C.OCOO 

P2#  0. 1680C000E-03 
P1**0. 36283040EC00 
PO#  0.  12266S63EE02 
r  3TK  0. 

12#  0. 

Tl#  C. 

TC#  0. S0879000EC03 
GC  TC  SC 

12  IF*H-?7COOO.nl3,  H,  |S 

13  X#*H-l7S00C.o/10000. 

P  3 #  -O.SOS7171 IE-03 
P 2#  -0. 71SS8221F-02 
PI#  -0.  3M69S39EC00 
PO#  0. I1723187EC02 
T3#  0. 12l90f 79ECC0 


ALT 
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LIST  OF  FORTRAN  PROGRAM 


SUPROLT INF 


T2#  -0.  13745684EG01 
TI«  -0  20357798EG02 
TOO  0. S075S1S2EC03 
GC  IL  5C 

14  CONTINUE 

15  IF*H-290000.ol6, 17,  18 

16  CONTINUE 

17  X0XH-27C000.O/100C0. 
P3#  0. 

P2*  0.26095006E-03 
PU-0.61251035EC00 
PO*  0*  7 1906  764EG01 
T3A  0. 

T2*  0. 

T  1  A  v. 

TOUT  0. 2982COOOE C03 
GO  TO  5 C 

IB  IFZH-35C00C. 019,20,21 
i  n  x*Xb-:> )  ;uo •: .  -i/ir.ooo . 

P  3  •  -0.  1224  30 1 4E-0 3 


al  r 
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LIST  OF  FORTRAN  PROGRAM 


SURROUT INC 


P 2#  0.  16998S82E-C1 

PI#  -0.  (S3252S19ECC0 
PC*  0.S9A791SIEC01 
T3#  -0.SS833219ECC0 
T2#  0.A1892744ECC1 

Tim  0.66193506EC00 
TO#  0.29A7S713EGC3 
GC  TO  SC 

20  CONTINUE 

21  X#*H-3S000C.o/lCCC0. 

P3#  -0. 263C9A06E-C2 
P2#  0. A  3  390852E-C 1 

PI#  -0.43767363CC00 
PO#  0.27SARS34EC01 

T 3#  -0. l0924334Ft00 
T2#  0. 8A07A3AIEC0O 

11#  0.  102S30S6EC03 

TO#  0. A03813A1ECC3 
SO  IFXH-29SO00. 022,22,23 

22  WM#28 • 9AA 


ALT 
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LIST  OF  FORJRAN  PROGRAM 


SUHROG  T I NE  AL1 


GC  UJ  SI 

2 3  IFth- 3SCPOC. a24,25,25 

24  WR#28.9AR8S8-tO.  1  1 7  1 4 744E-0  l n. X- tO . l42flSJ28F-02n*X*X 
GC  TC  SI 

25  MR #? 8. 84 892 7- tO.  26 9*2  1  llE-01u*X-t0.89303S08E-03u*X»X 
SI  PRriEXPFtttP3*XGP2o*Xt.Plu«X£POoo/l  00000. 

TEMIHttT  3<*XET2n«XCIlo»XGTO 
RG#tPR»V»M0/tlS4S.«TEMa 
RF  TORN 

SCS»SQRTFt4S.0416*PR/R0o 
FP# I . / l . 7406976ECC9-WM/RH 
VB»SCRTFt3.666M66rt.06*TLPii/WR 
CF#VB/FP 

GRA»  I.  3994IR2EUA/tHt2.08SSS3tCC7a*«2 

RETORN 

END 
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Line  -  0 
I  Page  =  1 


Print  106 


ENTER  FROM  PREVIOUS  PAOI- 
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ENTER  FROM  PREVIOUS  PAGE 


YR  =  RER  (COSF(ALA)) 
(SINF(ALO)) 


ENTER  FROM  PREVIOUS  PAGE 


X  =  RV  (COSF(PHI)) 
(COSF(THETA) 


ENTER  FROM  PREVIOUS  PAGE 


Enter  No.  1 
F rom  Page 
or  58 


ENTER  FROM  PREVIOUS  PAGE 


ENTER  FROM  PREVIOUS  PAGE 


ENTER  FROM  PREVIOUS  PAGE 


AK2  =  (AK1 2  +  2  (AK22)  + 
2(AK32)  +  AK42)  /  6 


ENTER  FROM  PREVIOUS  PAGE 


AK6  =  (AK1  6  +  2  (AK26)  + 
2(AKJ6)  +  AK46)/6 


ENTER  FROM  PREVIOUS  PACE 


ENTER  FROM  PREVIOUS  PAGE 


ENTER  FROM  PREVIOUS  PAGE 


ENTER  FROM  PREVIOUS  PAGE 


Next 

Page 
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ENTER  FROM  PREVIOUS  PAGE 


(Next  Page) 


SUBROUTINE  COOD 


RETURN 


END 


l 


SUBROUTINE  QUAD 


60 


SUBROUTINE  RAE 


SUBROUTINE  RAE 
(A,  B,  C,  XI,  Yl,  7.1) 

- , -  J 


_ i _ 

A  =  n/"(X'T)2~  I  (Yl)?  +  (/ 1)2 


A1  AN  (/.I) _ 

J(X1)<  7  (Y  1  +  ( /  1 ) 3 

. 01745 


AT  AN 


B  = 


XI 

Yl 


.01745 


ENTER  1  FROM 
PREVIOUS  PAGE 


>  0 


=  0 


=  0 


w 

B  =  90 


B  =  u 

*7 - 

iO  V  B 

\±7 

B  =  1F0  4  B 

RETURN 

END 


ENTER L  FROM 
PREVIOUS  PAGE 


<0 


. w 

B  =  180  +  B 


SUBROUTINE  COODI 


I 
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SUHHOUTINK  FU !  i  i 


65 


